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Abstract 

The evolution of the electromagnetic coupling, a, in the momentum-transfer 
range 1800 GeV 2 < -Q 2 < 21600 GeV 2 is studied with about 40 000 Bhabha- 
scattering events collected with the L3 detector at LEP at centre-of-mass energies 
y/s = 189 — 209 GeV. The running of a is parametrised as: 

(n 2 \ — a ° 

~ l-CAa{Q*)' 

where «o = ot(Q 2 = 0) is the fine-structure constant and C = 1 corresponds to the 
evolution expected in QED. A fit to the differential cross section of the e + e~ — > e + e~ 
process for scattering angles in the range | cos#| < 0.9 excludes the hypothesis of a 
constant value of a, C — 0, and validates the QED prediction with the result: 

C = 1.05 ±0.07 ±0.14, 

where the first uncertainty is statistical and the second systematic. 
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1 Introduction 



A fundamental consequence of quantum field theory is that the value of the electromagnetic 
coupling, a, depends on, or runs with, the squared momentum transfer, Q 2 . This phenomenon 
is due to higher momentum-transfers probing virtual-loop corrections to the photon propagator. 
This process of vacuum polarisation is sketched in Figure ^ In QED, the dependence of a on 
Q 2 is described as [2]: 

where the fine-structure constant, «o = a(Q 2 = 0), is a fundamental quantity of Physics. 
It is measured with high accuracy in solid-state processes and via the study of the anomalous 
magnetic moment of the electron to be l/«o = 137.03599911±0. 00000046 [1]. The contributions 
to Aa(Q 2 ) from lepton loops are precisely predicted [3], while those from quark loops are 
difficult to calculate due to non-perturbative QCD effects. They are estimated using dispersion- 
integral techniques [4] and information from the e + e~ — > hadrons cross section. At the scale of 
the Z-boson mass, recent calculations yield a _1 (m|) = 128.936 ±0.046 [5]. Similar results, with 
smaller uncertainty, are found by other evaluations using stronger theoretical assumptions. For 
example, Reference 6 obtains a _1 (m|) = 128.962 ± 0.016. 

The running of a was studied at e + e~ colliders both in the time-like region, Q 2 > 0, and the 
space-like region, Q 2 < 0. The first measurement in the time-like region was performed by the 
TOPAZ Collaboration at TRISTAN for Q 2 = 3338 GeV 2 by comparing the cross sections of the 
e + e~ — > e + e~ and e + e~ — > e + e~fi + [i~ processes [7]. The OPAL Collaboration at LEP exploited 
the different sensitivity to a(Q 2 ) of the cross sections of the e + e~ — > e + e~ — > t + t~ and 

e + e~ — > qq processes above the Z resonance to determine a(37236 GeV 2 ) [8]. Information on 
a(m|) is also extracted from the couplings of the Z boson to fermion pairs [9]. 

Bhabha scattering at e + e~ colliders, e + e~ — > e + e~, gives access to the running of a in the 
space-like region. In addition, like other processes dominated by t-channel photon exchange, it 
has little dependence on weak corrections. The four-momentum transfer in Bhabha scattering 
depends on s and on the scattering angle, 6: Q 2 = t ~ — s(l — cos6>)/2 < 0. Small-angle and 
large-angle Bhabha scattering allow to probe the running of a in different Q 2 ranges. 

LEP detectors were equipped with luminosity monitors, high-precision calorimeters located 
close to the beam pipe and designed to measure small-angle Bhabha scattering in order to 
determine the integrated luminosity collected by the experiments. The L3 collaboration first 
established the running of a in the range 2.10 GeV 2 < —Q 2 < 6.25 GeV 2 [10] by comparing 
event counts in different regions of its luminosity monitor. More recently, the OPAL Collabo- 
ration studied the similar range 1.81 GeV 2 < -Q 2 < 6.07 GeV 2 [11]. 

The running of a in large-angle Bhabha scattering was first investigated by the VENUS 
Collaboration at TRISTAN in the range 100 GeV 2 < -Q 2 < 2916 GeV 2 [12]. Later, the L3 
Collaboration studied the same process at y/s = 189 GeV for scattering angles 0.81 < | cos^l < 
0.94, probing the range 12.25 GeV 2 < -Q 2 < 3434 GeV 2 [10]. 

This Letter investigates the running of a by studying the differential cross section for Bhabha 
scattering at LEP at y/s = 189 — 209 GeV for scattering angles such that | cos#| < 0.9. Less 
than 1% of the events scatter backwards, cos 6 < 0, and this analysis effectively probes the 
region 1800 GeV 2 < — Q 2 < 21600 GeV 2 , extending and complementing previous space-like 
studies. 
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2 Analysis Strategy 



In the following, the running of a is described by a free parameter, C, defined according to: 



where the parametrisation of Reference 5 is used for the term Aa(Q 2 ). A value of C consis- 
tent with C = 1 would indicate that data follow the behaviour predicted by QED, while the 
hypothesis a = «o, with no dependence on Q 2 , corresponds to C — 0. 

The value of C is derived by a study of the measured differential cross section of the 
e + e~ — > e + e~ process, da /d cos 9. This quantity depends on C through the measured integrated 
luminosity, C{C), which is calculated from the expected cross section of the e + e~ — > e + e~ 
process for small scattering angles. The measurements used in the following are obtained under 
the Standard Model hypothesis, C = 1, as: 

MX) N(cos9) 1 

dcos# Acos# £(l)e(cos0)' 1 ' 

where N(cos9) is the number of events observed in a given cos 9 range, of width Acos#, with 
average acceptance e(cos9). The measured integrated luminosity depends on C as: 

C{C) = *l(C)s l (C) = C{1) a L L {C)e L {Cy (4) 

where Nl is the number of events observed in the fiducial volume of the luminosity monitor, 
&l(C) is the corresponding e + e~ — > e + e~ cross section for a given value of C and £l(C) is the 
detector acceptance. This acceptance may depend on C due to the combined effect of small 
angular anisotropics of detector efficiencies and the dependence of the predicted differential 
cross section on C. These changes in the acceptance are found to have negligible impact on the 
results presented below. 

The value of the parameter C is extracted by comparing the measured differential cross 
section to the theoretical prediction as a function of C, da th (C)/dcos9, derived as: 

da th (C) = da th (l) £(1) 

dcostf = dcostf £{Cy ^ ' 

where da th (1) / d cos 9 is the Standard Model prediction, discussed in Reference 13. The value 
of £(1) is derived by using the BHLUMI Monte Carlo program [14]. The dependence of 
da th (C)/dcos9 and C(C) on C is implemented by means of the BHWIDE Monte Carlo pro- 
gram [15]. The differential cross section is factorised as: 

"3— 4^ = —3 ^Frad COS 9) , 6 

d cos 9 d cos 9 

where da Born (C)/dcos9 is the tree-level differential cross section, which has a simple analytical 
form. The term F rad (cos 9) parametrises initial-state and final-state radiation effects, domi- 
nated by real-photon emission, as implemented in BHWIDE. It is verified that F Ta d (cos 9) has 
a negligible dependence on the spread of y/s considered in this analysis and, most important, 
on C. 
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3 Cross Section Measurement 



The data were collected at LEP by the L3 detector [16,17] in the years from 1998 through 2000. 
They correspond to an integrated luminosity of 607.4 pb" 1 and are grouped in eight intervals 
of y/s with the average values and corresponding integrated luminosities listed in Table [TJ 

Events from the e + e~ — > e + e~ process are selected as described in Reference 18. Electrons 
and positrons are identified as clusters in the BGO electromagnetic calorimeter, matched with 
tracks in the central tracker. In the barrel region of the detector, | cos^l < 0.72, the energy of 
the most energetic cluster must satisfy E\ > 0.25y/s, while the energy of the other cluster must 
satisfy E 2 > 20 GeV. In the endcap region, 0.81 < | cos#| < 0.98, these criteria are relaxed to 
Ei > 0.2y/s and > 10 GeV. Events with clusters in the transition region between the barrel 
and endcap regions, 0.72 < |cos#| < 0.81, instrumented with a lead and scintillating-fiber 
calorimeter [17], are rejected. To suppress contributions from events with high-energy initial- 
state radiation, the complement to 180° of the angle between the two clusters, the acollinearity, 
(, is required to be less than 25°. The number of events observed at different values of yfs is 
shown in Table Q together with the Monte Carlo expectations for signal and background. 

The e + e~ — > e + e~ process is simulated with the BHWIDE Monte Carlo generator as- 
suming (7=1. Background processes are described with the following Monte Carlo genera- 
tors: KORALZ [19] for e+e" -> r+r", KORALW [20] for e+e" -> W + W~, PYTHIA [21] for 
e+e" -> Ze+e", DIAG36 [22] for e+e" -> e+e'e+e", GGG [23] for e+e~ -> 777 and TEEGG [24] 
for e + e~ — > e + e~7 events where one fermion is scattered into the beam pipe and the photon 
is in the detector. The L3 detector response is simulated using the GEANT package [25], 
which describes effects of energy loss, multiple scattering and showering in the detector. Time- 
dependent detector inefficiencies, as monitored during the data-taking period, are included in 
the simulation. 

Systematic effects, such as charge confusion, are reduced by folding the differential cross 
section into d<r/d| cos#|, which is defined as: 

da dcr 1 da , , . 

j a \c.nnf)~>n' \ 1 



d| cos 6\ dcos6> lcose<0 dcos# lcose>0 
This differential cross section is measured in the fiducial volume defined by: 

12° < 6 C - )C+ < 168° (8) 
I cos0| < 0.9 (9) 
C < 25° (10) 

where 9 e - and 9 e + are the polar angles of the electron and the positron, respectively. The value 
of cos# is derived as: 

sin 16*0+ — O 



cosg= r ' e . ' . (11) 

sin u e - + sin u e + 

Ten intervals of | cos#| are considered for each of the eight values of y/s, for a total of 80 
independent measurements. Tableland Figure El present the measurements of dcx/d| cos#| and 
the Standard Model expectations. The larger uncertainties in the interval 0.72 — 0.81 are due 
to the transition region between the barrel and the endcap regions. 
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4 Results 



Figures H3 and 0] compare the combined differential cross section at the average centre-of-mass 
energy (a/5) = 198 GeV with the Standard Model prediction, corresponding to C = 1, and 
with a constant value of a, corresponding to C = 0. The data favour the hypothesis C — 1 
over the hypothesis C = 0, as also presented in Table El 

The value of C is extracted by comparing the 80 measurements of dtx/d| cos#| with the 
theoretical expectations da th (C)/dcos6 in a \ 2 fit with the result: 

C = 1.06 ±0.07, 

where the quoted uncertainty is statistical only. Several sources of systematic uncertainties are 
then considered. 

• The theoretical expectations for da th (l)/dcos8 have an uncertainty which varies from 
0.5% in the endcap region to 1.5% in the barrel region [13, 15]. 

• The measurements of der/d| cos(9| are affected by a systematic uncertainty, dominated by 
the event-selection procedure, which varies between 1% and 10%, as listed in TableEl [18]. 

• An uncertainty between 0.2% and 1.5% is assigned to F rad [cos 6), as a function of cos^, in 
order to account for possible higher-order effects not included in the BHWIDE parametri- 
sation. 

• Migration effects among the different cos 9 bins are found to be negligible due to the large 
bin size and the good detector resolution. 

Systematic uncertainties are conservatively treated as fully correlated and the fit is repeated 
including both statistical and systematic uncertainties with the result: 

C = 1.05 ±0.07 ±0.14, 

where the first uncertainty is statistical and the second systematic. A breakdown of the sys- 
tematic uncertainty is presented in Table EJ This result is in agreement with the Standard 
Model expectation, C = 1. The quality of the fit is satisfactory, with a x 2 of 91.9 for 79 degrees 
of freedom, corresponding to a confidence level of 17%. The hypothesis of a value of a which 
does not depend on Q 2 , C = 0, is totally excluded with a \ 2 °f 316 for 80 degrees of freedom, 
corresponding to a a confidence level of 10~ 29 . 

5 Discussion 

The result presented above establishes the evolution of the electromagnetic coupling with — Q 2 
in the range 1800 GeV 2 < —Q 2 < 21600 GeV 2 . This finding extends and complements studies 
based on small-angle Bhabha scattering by the L3 [10] and OPAL [11] Collaborations, which 
studied the regions 2.10 GeV 2 < -Q 2 < 6.25 GeV 2 and 1.81 GeV 2 < -Q 2 < 6.07 GeV 2 , 
respectively. The advantage of large-angle Bhabha scattering, investigated in this Letter, is to 
probe large values of — Q 2 , while studies of small-angle Bhabha scattering at lower values of — Q 2 
benefit from a larger cross section and thus statistical accuracy. The experimental systematic 
uncertainties of measurements in the two — Q 2 regions are implicitly different. At large — Q 2 , 
they are dominated by the selection of Bhabha events in the large-angle calorimeters, while at 
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low — Q 2 they mostly arise from the event reconstruction in the luminosity monitors and from 
effects of the material traversed by electrons and positrons before their detection. Both studies, 
at large and low — Q 2 , are affected by theoretical uncertainties on the differential cross section 
of Bhabha scattering, although in different angular regions 

Figures 03 and El present the evolution of the electromagnetic coupling with — Q 2 . A band 
for 1800 GeV 2 < — Q 2 < 21600 GeV 2 shows the 68% confidence level result from this analysis. 
It is derived by inserting the measured value of C with its errors in Equation (J2J) together with 
the QED predictions of Reference 5. The results from previous L3 data for Bhabha scattering 
at 2.10 GeV 2 < -Q 2 < 6.25 GeV 2 and 12.25 GeV 2 < -Q 2 < 3434 GeV 2 [10] are also shown. 
These two measurements are not absolute measurements of the electromagnetic coupling but 
differences between the values of ot(Q 2 ) at the extreme of the Q 2 ranges [10]: 

a~ 1 (-2.10 GeV 2 ) - « _1 (-6.25 GeV 2 ) = 0.78 ± 0.26 (12) 
a- 1 (-12.25 GeV 2 ) - oT 1 (-3434 GeV 2 ) = 3.80 ± 1.29. (13) 

The results in FigureElare obtained by fixing the values of a(— 2.10 GeV 2 ) and a(— 12.25 GeV 2 ) 
to the QED predictions of Reference 5 in order to extract the values of a(— 6.25 GeV 2 ) and 
a(— 3434 GeV 2 ) from Equations |T2|) and (|T3|) . The results shown in Figure El are obtained by 
first determining the values of a(— 2.10 GeV 2 ) and a(— 12.25 GeV 2 ) from the measured value of 
C and from Equation (j2J) and then extracting the values of a(— 6.25 GeV 2 ) and a(— 3434 GeV 2 ) 
from Equations ()12|) and (fT3*|) . This procedure relies on the assumption that the measured value 
of C also describes the running of the electromagnetic coupling for lower values of — Q 2 . Both 
figures provide an impressive evidence of the running of the electromagnetic coupling in the 
energy range accessible at LEP. 
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(VS) (GeV) 


^(pb- 1 ) 


N D 


Nmc 


N s 


N B 


188.6 


156.4 


11561 


11559 


11288 


271 


191.6 


29.7 


1976 


1953 


1905 


48 


195.6 


83.7 


5677 


5673 


5539 


134 


199.5 


83.5 


5382 


5338 


5201 


137 


201.8 


39.1 


2379 


2417 


2355 


62 


205.2 


75.9 


4259 


4165 


4063 


102 


206.7 


130.4 


7388 


7512 


7339 


173 


208.2 


8.7 


441 


484 


473 


11 


198.0 


607.4 


39063 


39101 


38163 


938 



Table 1: Luminosity-averaged centre-of-mass energies, (y/s), and corresponding integrated lu- 
minosities, £, used in the analysis. The y/s spread in each point is of the order of 1 GeV. The 
numbers of observed events, Nu, are given together with the total Monte Carlo expectations, 
NmC) an d their breakdown into signal, Ns, and background, Nb, events. The last row lists the 
average centre-of-mass energy, the total integrated luminosity and the total numbers of events. 
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dcr/d | cos#| (pb) 



(y/s) = 188.6 GeV 



(y/s) = 191.6 GeV 



(^) = 195.6 GeV 



(y/s) = 199.5 GeV 



cos 9 



(| cos 9 



Meas. 



Exp. 



Meas. 



Exp. 



Meas. 



Exp. 



Meas. 



Exp. 



0.00 
0.09 
0.18 
0.27 
0.36 
0.45 
0.54 
0.63 
0.72 
0.81 



0.09 
0.18 
0.27 
0.36 
0.45 
0.54 
0.63 
0.72 
0.81 
0.90 



0.052 
0.138 
0.227 
0.317 
0.407 
0.497 
0.588 
0.678 
0.770 
0.862 



12.8 ± 

10.9 ± 
14.0 ± 
16.2 ± 
25.0 ± 
35.0 ± 
57.9 ± 

109.8 ± 
227.1 ± 
735.4 ± 



0.9 ± 
0.8 ± 
1.0 ± 
1.0 ± 

1.3 ± 
1.6 ± 

2.0 ± 

3.1 ± 
18.2 ± 10.8 

8.4 ± 6.3 



0.2 
0.1 
0.2 
0.1 
0.2 
0.3 
1.1 
2.6 



10.4 
11.3 
13.4 
17.1 
23.7 
35.3 
57.7 
105.8 
232.2 
735.9 



9.3 ± 1.9 ± 0.2 

10.2 ± 2.0 ± 0.2 

11.5 ± 2.1 ± 0.2 

14.3 ± 2.4 ± 0.2 

23.6 ± 3.0 ± 0.3 
30.9 ± 3.5 ± 0.3 
61.2 ± 5.0 ± 1.2 

109.4 ± 7.4 ± 2.9 

196.4 ± 39.3 ± 16.3 

720.4 ± 19.8 ± 7.3 



10.1 
11.0 
13.0 
16.6 
22.9 
34.2 
55.9 
102.6 
225.1 
713.5 



8.6 ± 1.0 ± 0.1 

10.5 ± 1.2 ± 0.1 
12.0 ± 1.2 ± 0.2 
18.0 ± 1.5 ± 0.2 

20.6 ± 1.6 ± 0.2 
32.4 ± 2.1 ± 0.3 

51.2 ± 2.6 ± 1.0 

99.3 ± 4.0 ± 2.5 
211.2 ± 23.6 ± 14.0 
690.4 ± 11.2 ± 6.4 



9.6 
10.5 
12.4 
15.9 
21.9 
32.8 
53.6 
98.5 
216.2 
685.1 



10.3 ± 
9.8 ± 

12.2 ± 
14.7 ± 
20.0 ± 
28.0 ± 

49.3 ± 
98.9 ± 

231.3 ±26 
670.3 ± 11 



.1 ± 
.1 ± 
.3± 
.4± 
.6± 
.9± 
.6± 
.1 ± 



0.2 
0.1 
0.2 
0.2 
0.2 
0.3 
1.0 
2.5 



.9 ± 16.2 
.1± 6.2 



9.2 
10.0 
11.9 
15.2 
21.0 
31.5 
51.5 
94.6 
207.7 
658.4 



(y/s) = 201.8 GeV 



(y/s) = 205.2 GeV 



(y/s) = 206.7 GeV 



(y/s) = 208.2 GeV 



cos 9 



(\cos9\) 



Meas. 



Exp. 



Meas. 



Exp. 



Meas. 



Exp. 



Meas. 



Exp. 



0.00 
0.09 
0.18 
0.27 
0.36 
0.45 
0.54 
0.63 
0.72 
0.81 



0.09 
0.18 
0.27 
0.36 
0.45 
0.54 
0.63 
0.72 
0.81 
0.90 



0.052 
0.138 
0.227 
0.317 
0.407 
0.497 
0.588 
0.678 
0.770 
0.862 



11.0 ± 1.7 ± 0.1 

11.0 ± 1.7 ± 0.1 

11.9 ± 1.8 ± 0.2 

14.8 ± 2.1 ± 0.2 

21.2 ± 2.5 ± 0.2 

37.2 ± 3.3 ± 0.4 

55.5 ± 4.1 ± 1.0 

91.4 ± 5.7 ± 2.3 

243.7 ± 39.0 ± 18.7 

618.3 ± 15.8 ± 6.0 



9.0 
9.8 
11.6 
14.9 
20.6 
30.9 
50.5 
92.7 
203.7 
645.7 



8.7 ± 

12.9 ± 

12.3 ± 

16.1 ± 

20.0 ± 

31.7 ± 

48.0 ± 

93.3 ± 
252.2 ±29 
628.9 ± 11 



2± 
4± 
4± 
6± 
8± 
3± 
8± 
3± 



0.2 
0.2 
0.2 
0.1 
0.2 
0.4 
0.9 
2.3 



3 ± 15.5 
9± 5.7 



9.6 
11.4 
14.6 
20.2 
30.2 
49.5 
90.9 
199.7 
633.3 



9.0 ± 0.9 ± 0.1 

8.7 ± 0.9 ± 0.1 

10.4 ± 0.9 ± 0.2 

16.8 ± 1.2 ± 0.2 

23.1 ± 1.4 ± 0.3 

29.4 ± 1.6 ± 0.4 

44.5 ± 2.0 ± 0.8 
90.0 ± 3.2 ± 2.4 

170.0 ± 17.5 ± 14.5 

604.4 ± 8.7 ± 6.2 



8.6 
9.4 
11.2 
14.3 
19.8 
29.7 
48.5 
89.2 
195.9 
621.2 



3.9 ± 2.3 ± 0.1 

10.3 ± 3.7 ± 0.2 
5.4 ± 2.7 ± 0.1 

13.4 ± 4.2 ± 0.1 
23.4 ± 5.7 ± 0.3 
26.3 ± 6.0 ± 0.3 
37.6 ± 7.2 ± 0.7 
84.3 ± 12.0 ± 2.3 

280.3 ±88.6 ±24.0 

565.3 ± 33.0 ± 5.8 



8.5 
9.2 
10.9 
14.0 
19.4 
29.1 
47.6 
87.5 
192.2 
609.5 



Table 2: Measured, Meas., and expected, Exp., folded differential cross sections for the eight average centre-of-mass energies, (y/s), 
and the ten | cos#| intervals, with expected average values (| cos#|). The first uncertainty is statistical and the second systematic. 



<|COS0|> 



0.052 
0.138 
0.227 
0.317 
0.407 
0.497 
0.588 
0.678 
0.770 
0.862 



d I cos t 



(pb) 



9.93 ±0. 

10.25 ±0. 

11.99 ±0. 

15.95 ±0. 

22.15 ±0. 

31.65 ±0. 

51.15 ±0. 

98.7 ±1. 
211.6 ±9. 
666.9 ±4. 



42 ± 

43 ± 
47 ± 
54 ± 
64 ± 
77 ± 
99 ± 
5 ± 



0.15 
0.21 
0.14 
0.14 
0.25 
0.17 
0.26 
1.2 



± 13.9 
± 4.9 



d|cos0| yi JU ) 



9.7 
10.5 
12.4 
15.8 
21.7 
32.2 
52.3 
95.8 
210.2 
671.1 



do-"'(Q) 
d|cos6>| 



i (Pb) 



8.6 
9.4 
11.0 
14.2 
19.7 
29.5 
48.4 
89.1 
197.0 
634.2 



Table 3: Combined differential cross sections for the luminosity-averaged centre-of-mass energy 
(a/s) = 198 GeV, compared with the Standard Model expectations, dcr* /l (l)/d | cos#|, and the 
expectations for the case in which a does not change with Q 2 , da th (0)/d \ cos9\. The first 
uncertainties are statistical and the second systematic. 



Source of uncertainty 


AC 


Theoretical uncertainty 


0.11 


Experimental systematic 


0.08 


Frad 


0.05 


Bin migration 


< 0.01 


Total 


0.14 



Table 4: Sources of systematic uncertainty and their effect, AC, on the determination of the 
C parameter. 
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Figure 1: t-channel Feynman diagrams contributing to Bhabha scattering. Diagrams with 
virtual-fermion vacuum-polarisation insertions generate an electromagnetic coupling a(Q 2 ). 
The sum of all diagrams including zero, one, two or more vacuum-polarisation insertions is 
denoted by the diagram to the left with the double-wavy photon propagator. 
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Figure 2: Measured Bhabha differential cross- sect ions for the ten | cos#| intervals used in the 
study as a function of the centre-of-mass energy y/s. The Standard Model predictions are 
represented by the solid lines. 
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Figure 3: Measured differential cross-section as a function of | cos#|. Data at different centre-of- 
mass energies are combined at the luminosity- averaged centre-of-mass energy (\/s) = 198 GeV. 
The predictions in case of a running electromagnetic coupling and for a constant value a = «o 
are also shown. 
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Figure 4: Ratio between the measured Bhabha differential cross section as a function of | cos#| 
and the Standard Model expectations including a running electromagnetic coupling. Data at 
different centre-of-mass energies are combined at the luminosity-averaged centre-of-mass energy 
(a/s) = 198 GeV. The inner error bars represent statistical uncertainties and the full error bars 
the sum in quadrature of the statistical and systematic uncertainties. The predictions for a 
constant value of a = a are also shown. 
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Figure 5: Evolution of the electromagnetic coupling with Q 2 determined from the present 
measurement of C for 1800 GeV 2 < — Q 2 < 21600 GeV 2 , yellow band, and from previous 
data for Bhabha scattering at 2.10 GeV 2 < -Q 2 < 6.25 GeV 2 and 12.25 GeV 2 < -Q 2 < 
3434 GeV 2 [10]. The open symbols indicate the values of Q 2 where a{Q 2 ) was fixed to the QED 
predictions [5] in order to infer the values of a{Q 2 ) shown by the full symbols. These QED 
predictions are shown by the solid line. 
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Figure 6: Evolution of the electromagnetic coupling with Q 2 determined from the present 
measurement of C for 1800 GeV 2 < — Q 2 < 21600 GeV 2 , yellow band, and from previous 
data for Bhabha scattering at 2.10 GeV 2 < -Q 2 < 6.25 GeV 2 and 12.25 GeV 2 < -Q 2 < 
3434 GeV 2 [10], full symbols. The solid line represent the QED predictions [5]. 
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